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Abstract 
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1 Introduction 



In this note we present a leading evaluation of the coefhcient function of the shape function The latter is 

also called structure function of the heavy flavours. The coefficient function allows relating the shape function 
— computed with a non-perturbative technique such as lattice QCD — to distributions in semi-inclusive heavy- 
flavour decays. We consider in particular the processes 

B^Xs+i (1) 

and 

B^Xu + l + v (2) 

in the hard limit 

Q > A. (3) 

The quantity Q is the hard scale of these time-like processes: Q = 2£ , where £ is the hadronic energy of the 
state X and A is the QCD scale[]. For a i3-meson at rest, v = (1; 0, 0, 0) , and the jet X flying along the minus 
direction (— z axis), the shape function is defined as 

^ {k+) EE {B iv) \hl5 (fc+ - iD+) K\B [v)). (4) 

This represents the probability that the 6-quark in the B-meson has momenta 

Pb = ruBV + k' (5) 

with any transverse and minus components and with given plus component 

fc; = fc+. (6) 

The static field hy [x) is related to the Dirac field of the beauty quark h {x) by 

&(x) = e-*'"«''-"/i„(x)+o('— V (7) 



With the shape function, the decays (|l|) and (||) are related to their respective quark-level processes 

b^Xs + 1 (8) 

and 

h-^Xu + l + u, (9) 

where the 5-quark has the momentum (j^) with the distribution (^. It can be shown that the invariant mass m 
of the state X is related to the virtuality of the 6-quark by 

K - (10) 

The shape function is a non-perturbative distribution — analogous to parton distribution functions — and it 
describes the slice of the semi-inclusive region|^ in which 

~ AQ, 



^ For the rare decay (hi) , one can actually set Q = mg. 

^This region is also called threshold region, large-x region, radiation-inhibited region and Sudakov region. 
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i.e. 



|fc+|^A, or ^^1-^' (11) 



where 



(12) 

Because of ultraviolet divergences affecting its matrix elements, the shape function has a dependence on the 
ultraviolet cut-off or renormalization point /i : 

If {k+) = ip {k+; fi) , (13) 

and is related to a physical QCD distribution by means of a coefficient function by (cf. eq. (|77|)) 

cp{k+;Q) = j dk'^C{k+~k'^-Q,ti) ^{k'+■^J). (14) 

The QCD distribution does not depend on : 

j-cp{k+-Q) = Q, (15) 

while the shape function does not depend on Q. As a consequence, the coefhcient function depends on both Q 
and /i. 

The shape function can be computed with a non-perturbative technique or extracted from experimental data. 
If it is computed inside a field-theory model — such as lattice QCD — its expression will exhibit a /i-dependence 
that cancels against that of the coefficient function. If it is instead computed inside a phenomenological model — 
such as a quark model — the situation is less transparent. The /x-independence is not "automatic" and one has 
to specify the value of appropriate for the model. Some care is needed also in extracting the shape function 
from the experimental data, in order to avoid double counting of perturbative corrections. A factorization 
scheme must be defined and the coefficient functions for the various processes all have to be computed in the 
same scheme^. In particular, if a branching MonteCarlo is used for the analysis, the perturbative corrections 
generated by the program must be subtracted. 

The coefficient function is obtained by evaluating in leading approximation both and Lp and inserting their 
expression in eq. (|lj) . Since the coefficient function is expected to be a short-distance quantity, we compute 
the QCD distribution and the shape function in PT for an on-shell 6-quark {k' — 0). This expectation will be 
verified a posteriori. 



2 The QCD distribution 

The (perturbative) long-distance effects occurring in (|l|) and (^ can be factorized in the function 

f (z) = ^ (1 - z - 0) - A,as ) , (16) 



where 



M = ^ (17) 

TT 



^The situation is analogous to usual hard processes, where various factorization schemes for the parton distribution functions 
are defined: DIS, MS, etc. 
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and Cp — {N^ ~ l) / (2^c) = 4/3. The plus-distribution is defined as usual as 

P{z)^ = P{z) - Sil-z-0) f dyP{y). 

Jo 

The integrated or cumulative distribution is defined as 



F(z)= / dz'i{z'). 

J Z 

Inserting expression ( [l^ ) in this, one obtains the well-known double logarithm: 



Aias X 
——log [l- z). 



(18) 



(19) 



(20) 



The cumulative distribution satisfies the normalization condition F (0) = 1. Multiple soft-gluon emission expo- 
nentiates the one-loop distribution, so that 



F (z) ~ cxp 



Aias 2 



log' (1 - z) 



For further improvement, it is convenient to write the function f (z) in an "unintegrated" form, as 

i{z)=S{l-z) + A,as f- f' ^[S{l-z-et)-6{l-z)], 
Jo ^ Jo ^ 



where we have defined the unitary energy and angular variables 



E 

Q 



and t 



1 — cos ( 



(21) 



(22) 



(23) 



The quantity E is two times the energy of the soft gluon, E = 2Eg, and 9 is the emission angle. Leading 
logarithmic corrections are included replacing the bare coupling with the running coupling evaluated at the 
gluon transverse momentum squared 



with 



The QCD semi-inclusive form factor then reads: 

'■1 de dt 

T 



i{z)^5{l-z) + 
Performing the integrations, one obtains 



Aias (Q^eH) [S{l-z-€t)-5{l- z)] 



{{z)^S{l-z) + 



Po 



loglogC(l - z) 



1 - z 



log log ^(1 - zf 
1- z 



where /?q = (11/3 A^c — 2/3 n/) / (47r) and ^ is the square of the ratio of the hard scale to the QCD scale: 



A2' 



(24) 
(25) 

(26) 
(27) 

(28) 



Equation ( pTl) replaces the frozen-coupling result (|16[ ) and reduces to it in the limit ^ ^ oo with z fixed. The 
integrated distribution reads 



F (z) = 1 - [log ^ log log ^ - 2 log ^ (1 - z) log log ^ (1 - ^) + loge (1 - zf log log ^ (1 - zf 



(29) 
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As in the frozen-coupling case, higher-order corrections exponentiate the one-loop result, so that 

F (z) = cxp (-^ [log Clog log ^ - 2 log C (1 - z) log log e (1 - z) + log^ [l - zf log log ^ (1 ~ zf \. (30) 



Equation ( ^o[) replaces the frozen-coupling result . 

The above distributions are "physical" and therefore should be real for any value of z in the range 

0<z<l. (31) 
In practice, the distributions ( ^7|) and (|30| ) are real only if the range of z is restricted to 

z<l-^ or > KQ. (32) 

This restriction is absent in the frozen coupling case and originates from the infrared pole in the running 
coupling, which diverges when the transverse gluon momentum becomes as small as the QCD scale. The 
resummed distribution is therefore reliable as long as one stays away from the infrared pole, which implies the 
limitation 

m2>AQ. (33) 

This condition makes the kinematic region (^l]) unaccessible to the QCD distribution and forces the introduction 
of a non-perturbative component, namely the shape function. The physical origin of the restriction (^) is easily 
understood with the following qualitative considerations. The jet mass and the gluon transverse momentum 
are given by 

li^ Ed. (34) 

The smallest transverse momentum for fixed invariant mass is obtained for 6* ~ 1 and is 

TTl 

^±min ^ (35) 

In region this momentum is of the order of the QCD scale, where the coupling leaves the perturbative 
phase. The conclusion is that large-angle soft-gluon emission signals non-perturbative effects in the region (pi]). 

In general, resummation of multiple soft-gluon emission is performed in Mellin space. We therefore consider 
the moments of the form factor: 



1 

N , 



dz z'^ f (z) 

= 1 + Afjv. (36) 

Exponentiation of the "effective"^ one-gluon distribution takes place, so that 

f^^e"^^". (37) 

Using the large- A'^ approximation 

{1-yf -l^-e{y-l/n), (38) 



"^We call it "effective" because the insertion of the running coupling in the time-like region already includes some multiple 
gluon-emission effects. 
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where n = N/Nq with Nq = e '^^ ~ 0.56 and 7^ ~ 0.58 is the Euler constant, the moments of eq. jl^ ) or (22) 
read (frozen-couphng case, /3g — > 0) 



f„ = exp 



^i^s , 2 
r — log n 



In the running coupUng case, the MeUin transform of eq. (H) or is [|, 0, | 



f„ = exp 



? C C C 

log — log log — - 2 log - log log - + log ^ log log ^ 
n'' n n 



Note that f„ = 1 for n = 1. The distribution ( [40|) is usually written as 

f„ = exp gi iPoasl)] , 

where / = logn, as = as (Q) and 

gi (w) = - [(1 - 2w) log (1 -2w)-2{l- w) log (1 - w) 



(39) 



(40) 



(41) 



(42) 



Expanding this function to lowest order in w, one obtains gi {w) = — ^iw/(2/3q) + O {w'^) and recovers the fixed 
coupling result. 



As we clearly see, the inverse transform from f„ to F (z) is simply computed with the replacement [| 

1 



1 - z 



The limitation to the resummed perturbative result found before (eq. (|33[)) reads, in A^-space: 

Q 



n < 



A' 



(43) 



(44) 



3 The shape function 



Let us now consider the quantity in the effective theory related to f (z), namely the shape function ip (fc_|_) . For 
an on-shell 6-quark, the latter is given by 

^ (fc+) = 5 (fc+) + A,as ^^"^^t''^^ log ^ - A,as 5 (fc+) T ^ log (45) 

where /z = 2A5 and we defined 6 (ai; 02; • • ■ a„) = 9 (ai — 02) 9 (02 — 03) ■ ■ ■ 9 (a„_i — a„) . The regularization 
used imposes a cutoff on the spatial loop momenta and not on the energies: 



\ I \ < A5, —00 < Iq < +00; 
it is qualitatively similar to lattice regularization. Proceeding in a similar way, one can write 



(fc+) = 6 (fc+) + y^i«5 (Eh) [S (fc+ + i^t) - <5 (fc+)] 



We explicitly see that the hard scale Q does not appear in / (fc+), as it should. Since 
shape function in the "QCD variable" z reads 



/W = '5(] 



1 - z) + r - ^A,as (Q^eh) [S [1 - z - et) - S (1 - z)] , 
Ja ^ Jo t 



(46) 
(47) 

'Q (1 - z) , the 
(48) 
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where we have defined an adimensional shape function as 

The quantity 77 is the ratio of the UV cutoff of the effective theory to the hard scale, 



(49) 



(50) 



because the shape function is defined in a low-energy effective theory. To avoid substantial finite cut-off effects, 
it must also be assumed that 

M>A. (51) 
Taking the Mellin moments and exponentiating the one-loop distribution as in the QCD case, one obtain^ 

/w = e^/". (52) 

A straightforward computation gives 



log/„ 



A 



log — log log — - 2 log — log log h log 77^^ log log 77^^ 

n n 



(53) 



As in the case of the QCD distribution, /„ — 1 for n — 1. The A^-moments of the shape function are reliably 
computed in PT if one restricts n as in eq. (^4|). We therefore have here the same limitation of the QCD 
distribution, as expected on physical ground. 



4 The coefficient function 

We now introduce the coefficient function C (z) , relating the shape function to the QCD semi-inclusive form 
factor, as 

C{z) = 6{l- z) + AC{z), (54) 

where 

AC{z)=f{z)- f{z). (55) 
Inserting the expressions for the two distributions, we obtain 

AC (z) ^ r — r -A,as [Q^eh) [5 {1 - z - et) - 6 {1 - z)] . (56) 

Jri ^ Jo t 

The expression for the coefficient function is similar to the one for the QCD form factor, the only difference being 
a lower cut-off on the gluon energies in this case. Note that there is no angular cutoff. The QCD distribution 
and the shape function are related by a convolution in momentum space, 

i{z)^ f f dz'dz"d{z- z'z") C{z') f(z"). (57) 
Jo Jo 

Taking the moments on both sides, one diagonalizes the convolution, so that 

tn=Cnfn. (58) 



^Do not confuse these moments of the shape function /jv with the moments considered by other authors, J dk^k^Lp{kj^) 

fdz{i-zrf{z). 
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The usual exponentiation reads: 



Cn — e 



(59) 



The coefficient function can be computed from the difference in eq. ( pq ) or directly from the integral represen- 
tation in eq. (pq) in the following way. The integration over the emission angle t gives 



y 



max(?7, y) 



(60) 



where y = 1 — z. The integration over y is done using the previous approximation (|3^ ) and one finally obtains: 



logC„ 



+0(71-1/77) 



(1/77-71) 



log \ log log - 2 log - log log - + log ^ log log i 

n n n n 



+ 



(61) 



2 log — log log — — 2 log — log log — h log log log ^ — log 77^^ log log 77^^ 

71 71 77 77 



The QCD function f„ is equal to the coefficient of 9 — n) in the r.h.s of eq. (|6l|). Collecting the terms in 
common to the two 0-functions, eq. (61) can also be written as: 



logC„ = -^<j 6'(l/77-7i)log^loglog^ +log^loglog^-21og-loglog- + 



2/3o 



+9 (77-1/77) 



2 log — log log — - log 7;^^ log log 77^^ 
n n 



(62) 



Let us comment on the result, represented by eq. (^) or by eq. (£^). It is immediate to check that C„ = 1 when 
77 = 1 or when 77 = 1, as it should. Furthermore, C„ 1 in the limit Q — > 00 with 77 and fixed. 




Figure 1: Plot of the rcsummed coefficient function C„ (solid line) and the QCD semi- inclusive form factor f, 
(dashed line) for the choice of the parameters discussed in the text. 
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The main point is that the coefficient function is perturbatively computable as long as 

n < — -— . 

Note that this critical value of n is much larger than the one for the QCD form factor or the shape function 
given in eq. (^), because of eq. (|5l|). The equivalent limitation in momentum space for C (z) is 

m :» . (63) 

M 

The occurrence of non-perturbative effects outside the region indicated by eq. ( |63| ) can be understood in a simple 
way. The problem is to find the smallest gluon transverse momentum for a given invariant mass, eqs. (|^), with 
the additional constraint E > fx. The minimum occurs for the smallest gluon energy, E ^ fi, where one has 




and 



^ _Lmin 




(64) 



(65) 



in which conditions (^T|) and (|5l|) have been used. A comparison between eqs. (|3^) and ( |65| ) shows that trans- 
verse momenta are substantially larger in the coefficient function than in the QCD distribution; this allows a 
perturbative treatment of the former. This is because the infrared cut-off on the gluon energies in C (z) has the 
effect of lowering the largest emission angle, indirectly increasing the minimum gluon transverse momentum. 

In single logarithmic problems, the loop expressions for the coefficient functions usually have phase-space 
restrictions of the form 

> Ai > A. (66) 

Therefore, factorization with the shape function involves, as the relevant dynamical scale, a softer scale with 
respect to single logarithmic problems, because 

fi > ^/AJl. (67) 
Non-perturbative corrections to the factorized form ( |57| ) or (|5|) are expected to be of the size 

A2 



_Lmin 

which means of order 

A2 A 



(68) 



(69) 



One inverse power of the factorization scale ^ is involved. The corrections are therefore much larger than in 
single logarithmic problems, where the estimate (|68| ) translates to 

A2 

(70) 

involving two inverse powers of the factorization scale. 

The previous findings can be summarized by considering eq. (|5^) in the various regions of the moment index 
n, corresponding to different dynamical regimes: 
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1. In the regior 



1 « n « ^ 
A 



or equivalently 

Ag < m2 ^ (71) 

the QCD function, the shape function and the coefficient function are all reliably computed in perturbation 
theory. Actually, in this case there is no need to introduce the shape function and the resummed QCD 
form factor is all is needed: its splitting in a shape function and a coefficient function is irrelevant; 



2. In the "more exclusive" region 



or equivalently 



A A2 



AQ ^ » (72) 

the QCD function and the shape function become non-perturbative, while the coefficient function is still 
perturbatively computable. The idea is that one replaces the perturbative evaluation of the shape function 
with a non-perturbative one and defines the QCD distribution by means of relation (p8[); 

If one considers larger moments or, equivalently, smaller jet masses, 

n > -pj- or m < , (73) 

A-^ ^ 

also the coefficient function becomes non-perturbative and the shape function loses most of its meaning. 



The coefficient function C„ is plotted in fig. 1 together with the QCD function f„. For small values of n the 
two curves are very close to each other because typical transverse momenta are large and the lower cut-off ^ in 
Cn is ineffective. In the figure we have taken Q = 5.2 GeV, as = 0.24, /i = 2 GeV, A = 0.3 GeV and n/ = 3. 
For these value of the parameters, f„ becomes singular for n > 17, while C„ becomes singular for n > 115. The 
coefficient function monotonically decreases with n, starting from C„ = 1 at n = 1 and going down to C„ ~ 0.68 
for n = 17. This decay is produced by gluons of energy between the cutoff /i and the hard scale Q. The function 
f„ also monotonically decreases with n starting again from f„ = 1 at n = 1 and reaching f„ ~ 0.33 for n = 17. 
This faster decay is produced by gluons with any energy between the soft scale j/c+j and Q. 

The result (|6^) applies also to the coefficient function of the shape function defined in lattice regularization Q , 
after the identification is done 

- = c/i, (74) 
a 

where c is a constant of order 1 and a is the lattice spacing. The precise value of c is determined only in 
next-to-leading approximation. One has to perform a full order as computation of the shape function in lattice 
regularization and equate it to the expansion to order of the next-to-leading analogue of eq. (|4^). 

^The region n ~ 1 is trivial as there are no large infrared logarithms and fixed order perturbation theory is sufficient. 
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5 Conclusions 



In general, the QCD semi-inclusive form factor f (z) factorizes and resums large infrared double-logarithms 
which occur in the threshold region 

< (75) 

These terms are related to the emission of soft gluons with transverse momenta ranging from the infrared scale 
|fc+| = m? /Q up to the hard scale Q. As long as one keeps |fc+| ^ A, dynamics is controlled by perturbation 
theory. A phenomenologically relevant region is, however, |fc_|-| ~ A, in which substantial non-perturbative 
effects are encountered. We factorize them by means of the shape function. An additional, unphysical, scale /x 
is introduced in the problem, intermediate between the hard and the soft scale: 

|fc+| < A* < Q- (76) 
This means a splitting of the double logarithm in the QCD cumulative distribution of the form 



, Aias , 2 Q -, Mas , 2 Q 
1 — log — = 1 TT- log TT^ 



= 1- 



2 

Aias 



log' 



Ifc4 



Aiaslog — logT^ 

II fj, \k+\ 



Aios 2 ^J■ 



(77) 



The scale n has the role of an IR cut-off for the coefficient function and of an UV cut-off for the shape 
function. The first parenthesis in ( |7^ ) contains the coefficient function, which is short-distance-dominated and 
reliably computed in PT, while the second parenthesis contains a perturbative evaluation of the shape function. 
The latter computation is "unreliable" , so it is thrown away and replaced by a non-perturbative one. The 
factorization scale /x then separates what we compute in PT (above /i) from what we give up to compute in PT 
(below n). We may say that the coefficient function of the shape function is the "upper part" of a semi-inclusive 
QCD form factor. 

We have shown that the transverse gluon momenta entering the coefficient function have a lower bound 
given byQ 

li- > VW^ ~ » A, (78) 

where eq. ( [Ill ) has been used. As a consequence, non-perturbative corrections to the factorized form ( p7| ) or 
(H^) are expected to be of the size 

A , . 

involving one inverse power of the factorization scale /i. This situation is to be contrasted with single logarithmic 
problems, where corrections to factorization have typically a size 

A2 

(80) 

involving two inverse powers of the factorization scale. We may say that factorization in this double-logarithmic 
problem is still consistent, but it is presumably "less lucky" than in single-logarithmic cases. 

We conjecture that a similar factorization — into a perturbatively calculable coefficient function and a non- 
perturbative component — can be done for other semi-inclusive distributions, such as shape variables for very 

'^The appearance of a single logarithm of the soft scale |fc+| in the coefficient function, eq. ^n\l, should not be erroneously 
interpreted as the signal of a non-perturbative effect in C. 
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small values of the resolution parameters. An example is the thrust distribution in e+e annihilations |10(| for 
values of the jet masses in the region (0), for which 

(81) 

In other words, we believe that the idea of the shape function can be generalized to describe many semi-inclusive 
distributions in the region ([Tl|). 

To conclude, we have presented a leading evalutation of the resummed coefficient function of the shape 
function and we have proved that the latter is short-distance-dominated. 
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